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ABSTRACT 

We investigate the validity of the Cosmological Principle by mapping the cosmological parameters 
Hq and go through the celestial sphere. In our analysis, performed in a low-redshift regime to follow 
a model-independent approach, we use two compilations of type la Supernovae (SNe la), namely the 
Union2.1 and the JLA datasets. Eirstly, we show that the angular distributions for both SNe la 
datasets are statistically anisotropic at high confidence level (p-value < 0.0001), in particular the JLA 
sample. Then we find that the cosmic expansion and acceleration are mainly of dipolar type, with 
maximal anisotropic expansion [acceleration] pointing towards (/, h) ^ (326°, 12°) [(/, h) (174°, 27°)], 
and (/,6) (58°,—60°) [(/,0 (225°, 51°)] for the Union2.1 and JLA data, respectively. Secondly, 

we use a geometrical method to test the hypothesis that the non-uniformly distributed SNe la events 
could introduce anisotropic imprints on the cosmological expansion and acceleration. Eor the JLA 
compilation, we found significant correlations between the celestial distribution of data points and 
the directional studies of Hq and go, suggesting that these results can be attributed to the intrinsic 
anisotropy of the sample. In the case of the Union2.1 data, nonetheless, these correlations are less 
pronounced, and we verify that the dipole asymmetry found in the Hq analyses coincides with the 
well-known bulk-flow motion of our local group. Erom these analyses, we conclude that the directional 
asymmetry on the cosmological parameters maps are mainly either of local origin or due to celestial 
incompleteness of current SNe la samples. 

Subject headings: Cosmology: distance scale, cosmological parameters, Hubble diagram 


1. INTRODUCTION 

Spatial homogeneity and isotropy are large-scale prop¬ 
erties incorporated into the standard description of the 
Universe via the Cosmological Principle (CP) (see, e.g., 
Goodman, 1995; Maartens, 2011; Clarkson, 2012 for 
a discussion). According to current cosmological data 
(Eisenstein et al.^ 2005; Blake et al 2011; Suzuki et 
ai, 2012; Bennett et al, 2013; Ade et al, 2013; Betoule 
et al, 2014), the best picture nowadays is provided by 
the so-called ACDM model, a homogeneous and isotropic 
scenario which accounts well for the current accelerating 
phase of the Universe, most of the physical features in 
the Cosmic Microwave Background (CMB) spectrum and 
the observed large-scale structure. 

Erom the theoretical viewpoint, however, it is well 
known that there is still no satisfactory description of 
the mechanism behind cosmic acceleration which is a cru¬ 
cial aspect to the cosmological modelling. This in turn 
motivates the need to probe fundamental hypotheses in 
Cosmology which includes the validity of the CP. More¬ 
over, recent analyses have reported cosmological anoma¬ 
lies that may indicate deviations from an isotropic Uni¬ 
verse, such as large scale bulk-flow (Colin et al, 2011, 
Turnbull et al.^ 2012, Eeindt et al, 2013), alignment of 
low multipoles in the CMB power spectrum, CMB hemi¬ 
spherical asymmetry (Eriksen et al, 2004; Bernui et a/., 
2006), large scale alignment in the optical polarization 
QSO data (Hutsemekers et al, 2005; Hutsemekers et a/., 
2014) and spatial dependence of the value of the fine 
structure constant (Mariano & Perivolaropoulos, 2012). 
If an isotropic expansion is observationally confirmed, 
it would reinforce the search for a cosmological descrip¬ 
tion in the Eriedmann-Lemaitre-Robertson-Walker class. 


Otherwise, models with a more general spatial geometry 
should be seriously explored to describe the structure 
and evolution of the Universe. 

The idea of testing the CP with the Hubble diagram is 
recent. Previous efforts were done with several catalogs 
(Kolatt & Lahav, 2001; Schwarz & Weinhorst, 2007; An- 
toniou & Perivolaropoulos, 2010; Gumpta & Saini, 2010; 
Cai & Tuo, 2012; Kalus et a/., 2013; Zhao et al, 2013; 
Chang et al, 2014a, Gupta & Singh, 2014; Chang & 
Lin, 2015; Jimenez et a/., 2015), where the authors at¬ 
tempted to constrain the cosmological isotropy with the 
opposite hemispheres technique, or assuming anisotropic 
dark energy models (Blomqvist et al, 2010; Mariano 
& Perivolaropoulos, 2012; Cai et al, 2013; Yang et al.^ 
2014; Chang et al, 2014b; Wang et al, 2014). These ef¬ 
forts reported the validity of the CP within 2cr confidence 
level, with the dipolar directions of maximal anisotropy 
pointing in the direction of the cosmic anomalies we have 
mentioned earlier. 

The novelty of this work is a directional analysis of 
the cosmological parameters Hq and go in different sky 
patches by means of hemispherical comparison, including 
the investigation whether possible anisotropic signals in 
such analyses can be correlated with the non-uniform 
distribution of SNe events in the sky (according to the 
current catalogs). If a significant correlation is found, 
this would indicate that such signal could be attributed 
to observational biases rather than being of cosmological 
nature. We use the two latest compilations of SNe in 
order to perform these analyses, namely the Union2.1 
(Suzuki et al, 2012) and the JLA datasets (Betoule et al.^ 
2014) in a reasonably low-redshift regime, i.e., ^ < 0.20. 
This choice has been made as there are approximately 
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one third of the objects in both samples located within 
this range of redshift and also because of the advantage 
of carrying out model-independent analyses as pointed 
out by Seikel & Schwarz (2009) and Kalus et ai (2013). 
We use galactic coordinates throughout all our analyses. 


2 . SOME METHODS TO STUDY STATISTICAL 
ISOTROPY IN SNE DATA 

2 . 1 . Sigma-map 

In this section we describe our indicator and the pro¬ 
cedure to calculate their associated maps, which lead to 
quantify deviations from statistical isotropy in a given 
set of cosmic events with known positions on the celestial 
sphere (Bernui et a/., 2008; see also Bernui et a/., 2007). 
Our primary purpose is to illustrate the procedure for 
defining the discrete function a on the celestial sphere in 
order to generate their associated maps, called a—maps, 
which compared with statistically isotropic simulated 
maps give us a measure of deviations from isotropy in 
the data set. 

Let DU = D( 6 >j, 0^; 7 o) G 5^ be a spherical cap re¬ 
gion on the celestial sphere, of 70 degrees of aperture, 
with centre at the j-th pixel, j = l,...,Ycaps, where 
{0j^(j)j) are the angular coordinates of the center of the 
j-th pixel. Both, the number of spherical caps A^caps and 
the coordinates of their center {0j^(j)j) are defined using 
the HEALPix pixelization scheme (Gorski et a/., 2005). 
The spherical caps are such that their union completely 
covers the celestial sphere S ‘^. 

Let be the catalog of cosmic objects located in the j- 
th spherical cap DU. The 2 PACE of these objects (Pad- 
manabhan, 1993), denoted as Aj( 7 ^; 70 ), is the difference 
between the normalised frequency distribution and that 
expected from the number of pairs-of-objects with an¬ 
gular distances in the interval ( 7 ^ — 0.5(i, 7 ^ + 0.5J], i = 
1 , • • •, A^bins, where 7 * = (i - 0.5)(5 and 5 = 27o/A^bi„s is 
the bin width. The expected distribution is the aver¬ 
age of normalised frequency distributions obtained from 
a large number of simulated realisations of isotropically 
distributed objects in 5^, containing a similar number of 
objects as in the dataset in analysis. A positive (nega¬ 
tive) value of Aj indicates that objects with these an¬ 
gular separations are correlated (anti-correlated), while 
zero indicates no correlation. 

Let us define now the scalar function a : DU 3 ?+^ 
for j = I,..., Ycaps, which assigns to the j-cap, centered 
at ( 6 >j,(/)j), a real positive number (jj = (7{0j^(j)j) G 5R+. 
We define a measure a of the angular correlations in the 
j-cap as 

-j -^bins 

d = Y (1) 

’ bins • 

1 = 1 

To obtain a quantitative measure of the angular correla¬ 
tion signatures of the SNe sky map, we choose 70 = 90° 
and cover the celestial sphere with A^aps = Ahems = 768 
hemispheres, then calculate the set of values {dj, j = 
l,--,Vaps} using eq. O- Patching together the set 
{(Jj} in the celestial sphere according to a coloured scale 


(where, for instance, ^ blue^ ^maximal ^ 

we obtain a sigma-map. Einally, we quantify the angular 
correlation signatures of a given sigma-map by calculat¬ 
ing its angular power spectrum. Similar power spectra 
are constructed, for comparison, with isotropically dis¬ 
tributed samples of cosmic objects. 

2.2. Spherieal harmonies expansion 

Since the sigma-map assigns a real number value to 
each pixel in the celestial sphere, that is a = cr(^, 0 ), 
one can expand it in spherical harmonics: cr(^, 0 ) = 
Yj£,m^^rny£m{0,(t)) where the set of values {Q}, de- 

fined by Ci = (1/(2£ + 1)) is the angular 

power spectrum of the sigma-map. Because we are in¬ 
terested in the large-scale angular correlations, we shall 
concentrate on {Q, -^= 1 , 2 ,..., 10 }. 

2.3. The Huhhle- and q-maps 

This test involves the directional analyses of the Uni¬ 
verse expansion through the mapping of Hq and qo in 
the celestial sphere. In other words, we adopt the oppo¬ 
site hemisphere method to accomplish this. Each pair of 
these hemispheres is well defined by the HEALpix pix¬ 
elization scheme, such that we fit Hq and qo according 
to the maximal likelihood technique (MLT) in each one 
of them, and thus constructing the Hubble-maps and q- 
maps, respectively. The MLT consists in minimising the 
quantity 

, (2) 

i V / 

where the set (^^,/ii,cr^J contains the observational in¬ 
formation of the SNe data, i.e., redshift, distance mod¬ 
uli and associated uncertainty of the i-th object, respec¬ 
tively Q and /ith(^,p) is the distance modulus 

p) = 5 logio [Dl{z, p)] + 42.38 - 5 logio(h), (3) 

where h = Aq/IOO, Hq = lOOKm/s/Mpc, and Dl{z,P)) 
is the adimensional luminosity distance, whose argu¬ 
ments are the redshift z, in addition to the set of cos¬ 
mological parameters p which describe the underlying 
cosmological model. 

2.4. Statistieal signifieanee tests 

Once the Hubble- and g-maps have been constructed, 
it is of great relevance to test their statistical significance. 
Then, we follow two different approaches to accomplish 
this. In the first test one keeps the original set ( 2 :, /i, a^) 
of each SNe and shuffles their galactic coordinates, as 
performed, e.g., by Kalus et al. (2013). This proce¬ 
dure will be called the shuffle test throughout this paper, 
whose objective is to check the dependence between each 
(zi, /ii, with the SNe angular distribution pattern in 
the sky. Assuming that a SN event can occur in any 
patch of the sky independently of its redshift, it is not 
expected a dramatic change in the hemispherical asym¬ 
metries unless the cosmic expansion and/or acceleration 

^ is the error in the measurement of the SNe distance moduli, 
not to be confuse d wit h the discrete function cr^, j = 1, • • •, A^hems, 
defined in section [2.II to construct the sigma-map. 
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Figure 1. Left panel: The Union2.1 SNe distribution in the celestial sphere for 2 ; < 0.20 (230 data points). Right panel: The same catalog 
with the mask20 applied. 19 SNe have been excluded in this cut. 




Figure 2. Left panel: The Union2.1 SNe distribution in the celestial sphere for 2 ; < 0.20 (317 data points). Right panel: The same catalog 
with the mask20 applied. 15 SNe have been excluded in this cut. 


is really anisotropic, or the precision of the dataset is too 
limited. 

The second test also keeps the original (z, /i, <j^), how¬ 
ever, the SNe positions are randomly chose so the SNe 
are isotropically redistributed on the celestial sphere. 
This test will be referred as MC from now on. In this 
case, we are able to verify the constraining potential of 
our anisotropy diagnostics in the case of idealistically 
isotropic datasets. If any anisotropic signal persists, it 
must be accounted to a possible violation of the CP, or 
the limitation of the observational sample. 

Finally, we should remark that, although MC and shuf¬ 
fle can be carried out for the Hubble and q-maps, only 
the former is actually suitable for the sigma-map. Since 
the sigma-map only requires the angular position of the 
SNe events, the shuffle test is not convenient for these 
analyses as it maintains the original SNe distribution on 
the celestial sphere. Thus, as the sigma-map is designed 
to quantify the uniformity of angular distributions, only 
the MC realisations will be computed when the sigma- 
map test is performed. 

3. ANALYSES AND RESULTS 

3.1. The SNe eelestial maps 


Union2.1 sigma-map 


dipole amplitude 

{i,b) 

±0.0017 

(295.00°, -63.45°) 

JLA sigma-map 

dipole amplitude 

{l,b) 

±0.0020 

(195.00°, -81.22°) 


Table 1 

Amplitude and celestial position of the maximal sigma-map 
dipole contribution for the SNe Union2.1 (top), and the JLA 
datasets (bottom). Remember that this amplitude is given in 
arbitrary units, and that I G [0°,360°], b G [—90°, 90°]; for 
instance, the North (South) galactic pole has coordinates I 
arbitrary and b = 90° (b = —90°). In all the cases, the error in 
the angular estimates is ±3.66°. 


Considering only low-z SNe {z < 0.20) we plot the 
maps of the SNe distribution in the celestial sphere in 
figures [1] and [2] for the Union2.1 and JLA catalogs, re¬ 
spectively. In the top panel of both figures, it is possible 
to verify the more and less populated sky regions as well 
as a small number of SNe lying on the galactic plane due 
to the high thermal emission of dust in the Milky Way 
on that region. Hence, we also plot the SNe sky map 
with a galactic cut in the |6| < 20° plane, which can be 
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Figure 3. Left panel: Sigma-map of the low-z, mask20 Union2.1 catalog, showing their angular correlation signatures. The units are 
arbitrary. Right panel: The dipole term, i = 1, of this sigma-map, which indicates the direction in the sky where the mininal and maximal 
two-point angular correlations occur. 



Figure 4. Left panel: Sigma-map of the low-z, mask20 JLA catalog, showing their angular correlation signatures. The units are arbitrary. 
Right panel: The dipole term, ^ = 1, of this sigma-map, which indicates the direction in the sky where the mininal and maximal two-point 
angular correlations occur. 


visualised on the right panel of both figures. We name 
this galactic cut mask20^ and it will be used in our anal¬ 
yses in order to avoid any biases in our results due to 
incomplete coverage in this strip of the sky. 

It is worthy to remark that, even though the mask20 
cuts one third of the whole celestial sphere area, there is 
not so much data there due to the high optical extinction 
on this galactic region. The full-sky Union2.1 sample, 
limited to z < 0.20, contains 230 data points whereas the 
masked dataset has 211 events. It accounts less than 10% 
of the whole sample, therefore, it should not considerably 
affect the results of our analyses. For the JLA sample, 
the same procedure diminishes the number of SNe from 
317 to 302. 

3.2. The sigma-maps results 

The results for the sigma-map adopting the 2; < 0.20 
range and the mask20 galactic cut are presented in fig¬ 
ures [3] (Union2.1) and|J (JLA). The left panels of both 
figures show the original sigma-map, while the right pan¬ 
els feature their respective dipole-only contributions. Ta¬ 
ble 1 displays the amplitude of this dipole, in addition to 
the celestial position where it occurs in galactic coordi¬ 
nates. It is worthy remarking that the choice of N^ide = 8 


has been made for our plots so that the celestial sphere 
has 768 centres of hemispheres. 

From these analyses, we can readily note that both 
datasets are not perfectly isotropic, specially in the JLA 
sample, as the 2PACF yields significantly higher num¬ 
bers in some regions of the sky than others, pointing 
out a potential observational bias for anisotropic signal 
due to sky patches with under-sampled and over sampled 
SNe data. We measure the degree of anisotropic distri¬ 
bution of the SNe data by comparing its angular power 
spectrum with the mean of 500 sigma-maps power spec¬ 
tra produced by MC realisations, that is the case where 
the data points are isotropically distribuited in the sky. 
These results are shown in the figure [5] for both SNe cat¬ 
alogs. As the galactic cut encompasses one third of the 
whole celestial sphere, these simulated datasets contain 
33% additional events compared to the original sample, 
which would be the roughly expected number of observed 
SNe if the extinction around the galactic plane was much 
less intense. 

The goodness-of-fit test between the SNe and these 
simulated data gives us a measure of the consistency or 
deviation from the statistical isotropy hypothesis (that 
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Figure 5. Left panel: The angular power spectrum, {C^}, up to ^ = 10, of the sigma-map for the Union2.1 catalog. Right panel: the 
same for the JLA catalog. The crosses represent the values for the original datasets, while the red circles assign the average spectra from 
sigma-maps obtained from 500 MC realisations. Their error bars are estimated using the median absolute deviation of each coefficient of 
these spectra. The values, for 9 degrees of freedom, between real and simulated data are 5.6 x 10^ and 3.5 X 10^ for the Union2.1 and 
JLA samples, respectively, indicating that both SNe datasets are anisotropic at a high statistically significant level. 


Union2.1 Hubble-map 


dipole amplitude 

(i,b) 

±0.015 

(326.25°, 12.02°) 

JLA Hubble-map 

dipole amplitude 

{i,b) 

±0.025 

(58.00°, -60.43°) 


Table 2 

Amplitude and celestial position of the maximal Hubble-map 
dipole contribution for the SNe Union2.1 (top), and the JLA 
datasets (bottom). Remember that this amplitude is given in 
arbitrary units. In all the cases, the error in the angular estimates 
is ±3.66°. 


Uuiou2.I q-map 

dipole amplitude 

{i,b) 

±0.44 

(174.38°, 27.28°) 

JLA q-map 

dipole amplitude 

{i,b) 

±0.54 

(45.00°, -51.26°) 


Table 3 

Amplitude and celestial position of the maximal q-map dipole 
contribution for the SNe Union2.1 (top), and the JLA datasets 
(bottom). Remember that this amplitude is given in arbitrary 
units. In all the cases, the error in the angular estimates is 
±3.66°. 


is, the null hypothesis) of the SNe dataset. We obtain 
= 5.6 X 10^ and = 3.5 x 10^, respectively, for 9 
degrees of freedom, indicating that the SNe data from 
both catalogs are anisotropic at a high confidence level. 
Moreover, this multipole expansion analyses reveal that 
the dipole terms are much higher than the higher £ order, 
thus showing the significance of their contributions. 


3.3. The hubble and q-maps results 

As described in section 2, we construct our maps by as¬ 
sociating a pixel value with the best-fit of the parameter 
of interest in each one of the 768 hemispheres. Then, we 
are able to determine which sky patch presents the high¬ 
est and lowest expansion rate (in case of Hubble-map), 
the same for the cosmological acceleration {qq case). As 
our analyses are performed in the 2; < 0.20 range, we can 
use the cosmographic expansion to compute cosmological 
distances, which means that we do not need to assume 
any model for the dynamics of the Universe. Hence, the 
luminosity distance of the Eq. m is written as 

-Dl( 2 /,p) = 2 /+ y2 , y = (4) 

where qq = -ao/{aoH^). 

We have truncated our expansion series in the second- 
order term similarly to the procedure of Kalus et al. 
(2013), as the redshift range is short enough that the 
inclusion of third or higher order terms do not yield un¬ 
ambiguous results for Dp, independently of the value for 
go- It is also worthy noting that, as we have a two di¬ 
mensional parametric space (i7o, go), we marginalise over 
one of these parameters in order to map the other one 
on the celestial sphere. Therefore, marginalising over go 
corresponds to the Hubble-map case while marginalising 
over Hq gives us the g-map. 

The results for the Hubble-maps are shown in figures [6] 
and [7| for the Union2.1 and JLA samples, respectively, 
whereas the g-maps are displayed in figures [8] (Union2.1 
case) and [9] (JLA). We found the maximal asymmetry 
in the Hubble-map of SHq = iJomax “ ^Omin = 4.6 for 
the Union2.1 and SHq = 5.6 for the JLA case. This 
is larger than the measurement of asymmetry reported 
by Kalus et al (2013), i.e., SHq = 2.6, who, however, 
used a smaller SNe sample (183 objects) and fixed the 
value of the decelaration parameter at go = —0.601 in 
their analyses. The g-maps analyses yield Sqq = 2.56 
and Sqo = 1.62 for the Union2.1 and JLA SNe samples, 
respectively. 

The left panels of these figures present the dipole-only 
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Figure 6. Left panel: Hubble-map, in units of h = ifo/100, for the low-z, mask20 Union2.1 subset, showing the cosmological expansion 
in the celestial sphere. The lowest and highest values obtained for all 768 hemispheres are h = 0.683 and h = 0.729, respectively. Right 
panel: The dipole-only contribution term of the sigma-map, which indicates the direction in the sky where the mininal and maximal values 
of h occur. 



Figure 7. Left panel: Hubble-map, in units oi h = Hq/IOO, for the low-z, mask20 JLA subset, showing the cosmological expansion in 
the celestial sphere. The lowest and highest values obtained for all 768 hemispheres are h = 0.695 and h = 0.750, respectively. Right panel: 
The dipole-only contribution term of the sigma-map, which indicates the direction in the sky where the mininal and maximal values of h 
occur. 


contribution of each of these maps. These analyses al¬ 
low us to compute the preferred axes of our anisotropy 
diagnostics, that is Hq and qq. In the Hubble-map case, 
for instance, we obtained the highest expansion rate of 
the Universe towards (/,6) = (326.25°, 12.02°) for the 
Union2.1 catalog and (/,6) = (58.00°,—60.43°) for the 
JLA sample, as described in Table 2. For the q-map 
case, shown in Table 3, the preferred direction for the 
cosmological acceleration is (/,6) = (354.38°,—27.28°) 
(Union2.1) and (/,5) = (225.00°,—51.26°) (JLA). Even 
though there are some interception between both sam¬ 
ples, our Hq and qq analyses give rather different results 
for the amplitude of both anisotropy diagnostics as well 
as the directions in the sky where they occur. 

For completeness, we also compare our results with 
the preferred axes reported from different analyses (us¬ 
ing the Union SNe compilation) in Table 4. Finally, we 
also stress that this results seem to be consistent with 
the bulk-flow motion measurements reported by Turn- 
bull et al. (2012), who detected a bulk-flow of 249 ± 76 
km/s towards (/, h) = (319° ± 18°, 7° ± 14°) using nearby 
SNe, as well as Cohn et al. (2011), whose bulk-flow of 


(/,5) = (299°, 13°) was computed using the Union2 cat¬ 
alog at 0.015 < z < 0.2, the same redshift regime that 
we perform our analyses, and also Feindt et al (2013), 
who obtained the result of v^f = 292 ± 96km/s towards 
(/, b) = (290° ± 22°, 15° ± 18°) with the most nearby SNe 
(0.015 < 2: < 0.035) from this same sample. Thus, it is 
possible to identify the anisotropic signal of the Union2.1 
maps with this phenomenon. In the JLA case, however, 
we can not compare our results with previous reports 
since this sample contains many SNe events which are 
not present in the Union2.1 dataset. 

We also performed a simple test whose goal is to ver¬ 
ify whether the goodness of these fits coincide with the 
minimal and maximal values of Hq and qq that we have 
discussed earlier. Naturally, reduced i.e., xt = 
is adopted for these tests, where u represents the num¬ 
ber of degrees of freedom in each hemisphere. We note 
that the poorest fits for each parameter range from ^0.6 
to 1.3 in the Union2.1 maps and ^ 0.45 to 0.75 in the 
JLA case, where these values do not coincide with the 
minimal and maximal Hq or qq in any case. Thus, we 
conclude that the MLT we have carried out do not in- 
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-18.2' 


5 + 11 . 2 ° 


- 11 . 2 ° 
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(326.25' 

(58.00°, 

(354.38' 

(225.00' 


^ 12 . 02 °) 
-60.43°) 
^ 27.28°) 
^51.26°) 


Ref^ 

Antoniou & Perivolaropoulos, 2010 
Cai & Tuo, 2012 
Kalus et al, 2013 
Zhao et al.^ 2013 
Cai et a/., 2013 
Yang et a/., 2014 
Chang et al.^ 2014 
Wang et al.^ 2014 
this paper (hubble-map, Union2.1) 
this paper (hubble-map, JLA) 
this paper (q-map, Union2.1) 
this paper (q-map, JLA) 


Table 4 

The axes reported in the literature, and our results, along which 
the maximal anisotropy occurs. 


troduce a bias on the coslogical expansion (or accelera¬ 
tion) signals. Similar analyses have also been performed 
by Kalus et al. (2013), where no impact between their 
Hubble-maps and their y^-maps have been detected as 
well. 


3.4. The eorrelation analyses 

In Table 5 we display the correlation between each 
map, where the Pearson’s coefficient p have been adopted 
as a diagnostics for this purpose. For the Union2.1 
dataset, it is possible to note that the Hubble-map and 
g-map show a significant correlation: p —0.70, roughly 
meaning that in the sky regions where Hq is large q is 
low and negative, as expected by the definition of q. Yet, 
the most important cases are the correlation between 
Hubble-maps and g-maps versus the sigma-maps. De¬ 
spite the visual correlation between the celestial regions 
containing a large number of SNe (Fig. [1]) and the re¬ 
gions where Hq achieves its lowest values (Fig. |4]), hence 
indicating a significant observational bias in our results 
as a matter of fact there is a small correlation between 
the Hubble-map as compared with the sigma-map, with 
Pearson’s coefficient p +0.059. This means that over¬ 
sampled regions do not necessarily influence the Hubble 
diagram result. Similarly, there is weak anti-correlation, 
p —0.200 between the g-map and the sigma-map. 
These analyses let us to realise an overall assessment 
that the anisotropic SNe distribution does not imply an 
a priori influence in the anisotropic Hubble- and g-maps 
results. Finally, we have also verified if there is a correla¬ 
tion between Hubble- and g-maps and their correspond¬ 
ing y^-map, finding out that this correlation is small in 
both cases. This result, therefore, confirms the validity 
of our approach regarding the absence of bias from the 
statistical technique we have adopted to map the cosmo¬ 
logical parameters in the sky. 

Nevertheless, these analyses provide distinct results for 
the JLA compilation, which can also be checked on table 
5. There is a relevant correlation between Hubble and 
g-maps with the sigma-maps, evaluated on p = +0.65 
and p = —0.45, respectively, besides a stronger correla¬ 
tion occurs between the Hubble and g-maps: p = —0.91. 
These results indicate that the incomplete sky coverage 
of this sample indeed affects the anisotropy of the cosmic 


Union2.1 maps 

P 

Hubble-map vs. sigma-map 

+0.059 

q-map vs. sigma-map 

-0.200 

q-map vs. Hubble-map 

-0.702 

Hubble-map vs. X^Hubble-map 

+0.054 

q-map vs. X+-maD 

-0.391 


JLA maps p 


Hubble-map vs. sigma-map 
q-map vs. sigma-map 
q-map vs. Hubble-map 
Hubble-map vs. X^Hubble-map 
q-map vs. X^-map 


Union2.1 vs JLA maps p 

sigma-map Union2.1 vs. sigma-map JLA +0.614 

Hubble-map Union2.1 vs. Hubble-map JLA —0.188 
q-map Union2.1 vs. q-map JLA +0.261 


+0.651 

-0.446 

-0.915 

-0.130 

-0.346 


Table 5 

The correlation between pairs of maps. We have adopted the 
Pearson’s coefficient p to measure the correlation between diverse 
pairs of maps. 


SNe catalogs 

MCs (%) 

Shuffies (%) 

Union2.1 (z < 0.2) hubble-map 

3.0 

7.6 

JLA (z < 0.2) hubble-map 

0.6 

18.2 

Union2.1 {z < 0.2) q-map 

0.0 

0.0 

JLA {z < 0.2) q-map 

4.2 

42.8 


Table 6 

The percentage of Shuffie and MC realisations whose Ah or Ag'o 
are larger than that of the original datasets. 


expansion and acceleration. 

Finally, we have also computed the correlation between 
the maps obtained for each dataset. Although there is 
a suggestive correlation between the sigma-maps, this 
trend is not followed by the Hubble- and g-maps, which 
are just mildly correlated. This analysis reflects the re¬ 
sults we have discussed before, where we have shown that 
the axes of maximal anisotropic expansion (and acceler¬ 
ation) are very different for each catalog. 

3.5. The multipole analyses and Monte Carlo tests 

The figures [10] and m features the multipole expan¬ 
sion for the Hubble-maps (left panel) and g-maps (right 
panel) for the original Union2.1 and JLA datasets, re¬ 
spectively, compared to 500 realisations of MC and shuf¬ 
fle tests, where table 6 features the percentage of Shuffie 
and MC realisations whose SHq or Sqo are larger than 
that of the the original datasets. 

In the left ofHO] we can note that, for most the orig¬ 
inal data present larger qs than the simulated datasets 
except for the dipole coefficient, whose concordance is 
mildly larger than Icr uncertainty. On the other hand, 
the asymmetry of iJo, dHo = 0.046, of these simulations 
does not present significant disagreement with the result 
obtained from data, as 7.6% of the shuffled realisations 
yield 6Hq > 0.046 of the Union2.1, whereas this number 
decreases to 3.0% for the MC test. By the same token, 
18.2% of the shuffie test exceed the asymmetry on JLA 
Hubble-map (0.055), however this percentage drops to 
0.6% when performing the MC realisations. This result 


























- 2.4 


0.16 


- 0.44 


Figure 8. Left panel: The q-map for the low-z, mask20 Union2.1 subset, showing the cosmological expansion in the celestial sphere. The 
lowest and highest values obtained for all 768 hemispheres are qo = —2.40 and qo = +0.16, respectively. Right panel: The dipole-only 
contribution term of the g'-map, which indicates the direction in the sky where the minimal and maximal values of q^ occur. 



Figure 9. Left panel: The q-map for the low-z, mask20 JLA subset, showing the cosmological expansion in the celestial sphere. The 
lowest and highest values obtained for all 768 hemispheres are qo = —2.35 and qo = —0.73, respectively. Right panel: The dipole-only 
contribution term of the ^-map, which indicates the direction in the sky where the mininum and maximal values of qo occur. 


also shows the intrinsic non-uniformity of the angular 
SNe distribution of this sample, as the MC test yields 
smaller asymmetric signals than the original map when 
we redistribute the objects in the sky in uniform manner. 

In the right panels of both figures, for instance, we 
stress the comparison between the g-maps and the sim¬ 
ulated realisations of the Union2.1 and JLA data. For 
the first compilation, the multipole expansion of the q- 
maps clearly indicates a significant asymmetry, as most 
^s (except for the dipole term) are not reproduced by the 
simulations we have produced, such that this asymmetry 
is reflected on the 6qo computation, where none of the 
MC and shuffle realisations exceed the value obtained 
for the original sample, that is Sqo = 2.56. Although this 
large asymmetry can not be reproduced by these statisti¬ 
cal tests, we notice that this discrepancy is mostly due to 
the presence of an outlier fit for go, he., gomax = +0.160, 
and that no other hemisphere provide a positive fits for 
the deceleration parameter except for this. Such large 
range of values allowed for go can be attributed to a low 
constraining power on go from this SNe dataset. Hence, 
it remains inconclusive whether the cosmic acceleration 
is actually anisotropic or not albeit the low statistical 


significance of this result. 

Nevertheless, the JLA q-map presents a better agree¬ 
ment with these statistical tests, specially the shuffle test, 
as shown in the multipole expansion on the right plot 
of the figure [TU in addition to the fact that 42.8% of 
these shuffle realisations present 6qo > 1.62. However, 
the MC test, show a smaller concordance with the data, 
as only 4.2% of the simulations are able to reproduce 
such anisotropy on the q-map. As expected, the smaller 
asymmetry on idealistically isotropic JLA samples rein¬ 
forces the impact of the incomplete sky coverage of these 
catalog. 

It is also worth to discuss that the directions of the 
maximal and minimal (iLo,go), besides the sigma-map 
cases, are randomly oriented in both statistical tests. 
These positions are reproduced by less than 2% of the 
500 realisations in all cases with a tolerance of ^ ±11° 
in both latitude and longitude coordinates. Therefore, 
it is possible to conclude that such directions does not 
present statistical significances, even though they coin¬ 
cide with other CP probes, and show suggestive concor¬ 
dance with the reports from the literature which have 
adopted similar methodologies. 
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Figure 10. Left panel: The angular power spectrum, {C^}, up to ^ = 10, of the Hubble-map for the Union2.1 catalog. Right panel: the 
same for the q-map. The crosses represent the values for the original dataset, while the red (blue) circles (squares) assign the average 
spectra from Hubble-maps obtained from 500 MC (shuffle) realisations. Their respective error bars are estimated using the median absolute 
deviation of each coefficient of these spectra. 
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Figure 11. Left panel: The angular power spectrum, {C^}, up to ^ = 10, of the Hubble-map for the JLA catalog. Right panel: the same 
for the q-map. The crosses represent the values for the original dataset, while the red (blue) circles (squares) assign the average spectra from 
Hubble-maps obtained from 500 MC (shuffle) realisations. Their respective error bars are estimated using the median absolute deviation 
of each coefficient of these spectra. 


To summarise, we cannot conclude that the anisotrop- 
ical signals found on the Union2.1 catalog analyses are 
influenced by the angular correlations of SNe distribu¬ 
tion in the sky, together with the reasonably low statis¬ 
tical significance of these results, however, the proximity 
of the bulk-flow motion directions reported in the liter¬ 
ature raises the possibility of being the cause for such 
signals. The JLA analyses, on the other hand, clearly 
show that its anisotropic coverage of the sky indeed im¬ 
pact the mapping of cosmological parameters through 
opposite hemispheres. 

Therefore, it is not possible to determine any violation 
of the CP, in this redshift range, with the limitation of the 
current datasets. Nevertheless, next-generation surveys 
such as LSST (Abell et al, 2009) and Euclid (Amendola 
et a/., 2012) may solve this issue with the greater preci¬ 
sion, and much larger datasets, that they shall provide. 

4. CONCLUSIONS 

In this work we have discussed the issue of a possible 
anisotropic expansion of the Universe through a direc¬ 
tional analyses, i.e., investigating the Hubble diagrams 


of SNe data located in opposite hemispheres. For com¬ 
pleteness, our analyses also explore the possibility that 
the results are biased by the anisotropic distribution of 
the SNe data, because they are manifestly under-sampled 
or over-sampled in several sky patches. For this, we 
have firstly used a geometrical method to find possible 
anisotropic signatures due to the angular distribution of 
both SNe datasets, namely the Union2.1 and the JLA 
catalogs. Then we use this information to investigate a 
possible correlation with the outcomes of the directional 
studies of Hq and q^. 

In fact, according to our analyses, a correlation exists 
between the anisotropic distribution of the data ensam- 
bles and the anisotropic signatures found in our direc¬ 
tional studies of Hq and go- Specifically, the JLA data 
shows a large correlation that is statistically significant 
to suggest that the dipolar direction of the Hubble-map 
can be attributed to the anisotropically distributed SNe 
sample. In the case of the Union2.1 catalog, this corre¬ 
lation is small. However, we have verified that the di¬ 
rection of the dipole asymmetry found in the Hq analy¬ 
sis with this dataset correlates well with the well-known 
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direction of the bulk-flow motion, in fact we obtained 
(/, 6) = (326°, 12°) while the bulk-flow motion of our lo¬ 
cal group points towards (/,6) = (319°, 7°) (Turnbull et 
al, 2012). 

All these results indicate that the dipolar features of 
the directional Hubble diagrams (which is the domi¬ 
nant anisotropic signal) can be explained either by the 
anisotropic SNe dataset (in the case of the JLA sample) 
or by a systematic effect, that is, the direction of the 
bulk-flow motion of our local group in the case of the 
Union2.1 catalog. Thus, for the low-redshift regime of 
the SNe analysed here {z < 0.2), we have shown that 
the observed anisotropic expansion cannot guarantee a 
violation of the CP. Despite of this, one cannot discard 
the hypothesis that our Universe is indeed undergoing 
an anisotropic expansion, which is not observed in the 
analyses because either it could be of small magnitude 
or it may be hidden by the uncertainties introduced by 
the highly non-uniform distribution of current SNe data 
in the celestial sphere. 


We are grateful to Roy Maartens, Ribamar Reis and 
Ivan Soares Ferreira for helpful discussions. We acknowl¬ 
edge CNPq, FAPERJ and CAPES for financial support. 
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